Instabilities of the Hubbard chain in a magnetic field 
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We find and characterize the instabilities of the repulsive Hubbard chain in a magnetic field by 
studing all response functions at low frequency u and arbitrary momentum. The instabilities occur at 
momenta which are simple combinations of the (U = 0) a =f , J, Fermi points, ±kFo- For finite values 
of the on-site repulsion U the instabilities occur for single a electron adding or removing at momenta 
±&Fct, for transverse spin-density wave (SDW) at momenta ±2/cf (where 2ftF ~ fej?f + kpi), and for 
charge-density wave (CDW) and SDW at momenta ±2fcF| and ±2fcj?i. While at zero magnetic field 
removing or adding single electrons is dominant, the presence of that field brings about a dominance 
for the transverse ±2&f SDW over all the remaining instabilities for a large domain of U and density 
n values. We go beyond conformal-field theory and study divergences which occur at finite frequency 
in the one-electron Green function at half filling and in the transverse-spin response function in the 
fully-polarized ferromagnetic phase. 

PACS numbers: 72.15. Nj, 74.20. -z, 75.10. Lp, 67.40. Db 



I. INTRODUCTION 



• One-dimensional (ID) electronic quantum problems have peculiar instabilities which follow from their restrictive 
geometry. For instance, it is well known that the low-frequency charge and spin response functions of a gas of electrons 
in a ID lattice show logarithmic divergences at twice the Fermi wavevector. This is the so called Peierls instability 
which in real quasi-one-dimensional materials occurs through coupling of the electrons to suitable phonon modes 
giving rise to lattice distortions and to phenomena like charge-density waves (CDW's) or spin-density waves (SDW's). 

These instabilities also occur in the presence of a magnetic field and, in real materials there are electron - electron 
interactions which at low dimensions have a strong impact on the physical properties. In this paper we study the effects 
of the interplay between a magnetic field and many-electron interactions on the ID quantum-problem instabilities. We 
consider an integrable many-electron problem , the Hubbard chain || in a magnetic field [||-|§ • One of our goals 
is to use elsewhere the present results on the ID instabilities as a starting point for the study of the corresponding 
instabilities of a system of weakly coupled chains. Therefore, we include in our present ID study the one-particle 
Green function which involves excitations that change the electron numbers. In the coupled-chain system this ID 
function will provide important information on the relevance of interchain electron hopping. In Sec. VI we will return 
to the problem of the Hubbard coupled-chain system in a magnetic field. 

At zero onsite interaction, U = 0, the Hubbard chain reduces to a ID tight-binding electron model with first- 
neighbor transfer integral t. In this case it is easy to evaluate the low-frequency response functions Rex^{k,uj) 
(with $ being charge p, spin projection a z , and so on) and to characterize and detect the corresponding divergences. 
At zero magnetic field the charge function Rex p (±2fcj?, ui), the spin function Rex CTz (±2fcj?, u>), the spin-transverse 
function Rex CTi (±2fcp, w), the singlet superconductivity functions Rex ss (0, oj) and Rex ss (±2fci?, w), and the triplet a 
superconductivity function Rex tsa (0, cj) all diverge as — ln(u). If we apply a magnetic field we find that the H = 
instabilities are limiting cases of the finite-field situation - for < H < — [sin^n)] 2 (here H c = — [sm(^n)] 2 is the 
U = critical field for onset of fully-polarized ferromagnetism and /iq is the Bohr magneton) the charge functions 
Rex p (±2fc_F f, lu) and Rex p (±2fcFj., w), the spin functions Rex' 7 * (±2fefi, uo) and Rex CTz (±2fci?|, w), the spin-transverse 
function Rex a± (±2fep, lo), the singlet superconductivity functions Rex ss (±[fcF| — fc-FjJjW) and Rex ss (±2kp,oj), and 
the triplet a superconductivity function Rex tscr (0,w) all diverge as — ln(u>). Also, both for H — and H > the 
single-electron Green function diverges at ±kp a as 1/ui. 
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On the other hand, the study of all divergences of the low-frequency response functions Rex^fc, u>) is for the 
corresponding many-electron quantum system a problem of some complexity. The BA solution and conformal-field 
theory (CFT) can be combined to evaluate the asymptotics of the correlation functions of the Hubbard chain in a 
magnetic field [^|-^). However, the momentum-dependent expressions derived in Ref. j^] (for zero magnetic field see 
also Rcf. ||) and some of the quantities studied in Ref. M refer to the equal-time situation and we need here the 
low-frequency response functions. Moreover, Ref. 0] evaluated and presented expressions for the response functions 
which refer to particular relations between the values of the low-frequency lu and small- momentum (k — ho). Here 
fco are particular combinations of the U = electronic Fermi momenta [see Eq. (12) below]. We find that the low- 
frequency correlation functions Hex^(k,uj) can only diverge at the particular values k = k$. Unfortunately, Ref. Q] 
has not solved the present problem of the instabilities because its expression (5.7) corresponds to a very particular 
relation between u and (fc — fc ) which does not provide the relevant and general oj dependence at k — k . 

The problem was solved for the particular case of the charge $ = p and spin ■§ = a z response functions in Ref. [[l0| . 
However, that paper has not considered all correlation functions. Although the charge and spin excitations studied 
in Ref. are relevant and correspond to real instabilities, they are never dominant, as we show in this paper. 

Using the same methods as in Refs. [ [iolJll) , we find that all instabilities correspond at finite values of U to power- 
law divergences (and in some singular points of the parameter space to logarithmic divergences) in the correlation 
functions. The corresponding exponents are non-classic combinations of two-pseudoparticle phase shifts. These 
pseudoparticles are the quantum objects introduced in Refs. [|l2-16| which refer to the operator algebra jlT] 
which diagonalizes the quantum problem. Their zero-momentum forward-scattering interactions control the overlaps 
between the Hamiltonian eigenstates and the excitations associated with the correlation functions 10 llj]. It is the 
exotic character of such overlap which leads to the Luttinger-liquid pfj| non-classic exponents. 

While at zero magnetic field removing or adding single a electrons at k = ikpa- is dominant (in the weakly coupled- 
chain problem this is associated with inter-chain hopping) and, due to the spin SU(2) symmetry, the transverse ±2kp 
SDW and the ±2kp SDW and CDW have the same oj dependence, the presence of a magnetic field brings about a 
dominance for the transverse ±2/c_f SDW over the remaining instabilities for a large domain of U and n values. There 
are no divergences in the superconductivity and other response functions. 

We study the effects of the half-filled metal - insulator transition and fully-polarized ferromagnetic transition on 
the above correlation functions. Some low-energy power-law divergences are present in both sides of these transitions 
but there occur discontinuities in the corresponding non-classical critical exponents. Moreover, we go beyond CFT 
and study power-law divergences which occur for some correlation functions at finite energy for the half-filling and 
fully-polarized ferromagnetic phases. 

The paper is organized as follows: the Hamiltonian and the pseudoparticle operator basis are introduced in Sec. 
II. In Sec III we characterize the instabilities of the quantum problem. The comparative study of these instabilities 
is presented in Sec. IV. In Sec. V we consider the half-filling and fully-polarized ferromagnetic phases and the effects 
of the half-filling metal - insulator transition and fully-polarized ferromagnetic transitions on the quantum-problem 
instabilities. The coupled-chain problem in the presence of a magnetic field is discussed in Sec. VI. Finally, Sec. VII 
presents the concluding remarks. 



II. THE HUBBARD CHAIN AND THE PSEUDOPARTICLE OPERATOR BASIS 



We consider the Hubbard chain [|5|-|7| with a finite chemical potential /i and in the presence of a magnetic field H 
ft = - *S [S> c i+ 1 - CT + c ]+i^A + U &i,T c i.T - l/ZHct^j. - 1/2] + 2^ z + 2fi HS z , (1) 

j,o 3 

where ct a and Cj. a are the creation and annihilation operators, respectively, for electrons at the site j with spin 

projection a =|, \ and i) z — —h[N a — J2 a ^<r] an d S z = — | J2a a ^a- Here N a is the cr-electron number operator. We 
introduce the U — Fermi points kp a — ~Kn a (and kp = [kp-\ + /cpjJ/2 = im/2) where n a = N a /N a and n = N/N a , 
and N a is the number of lattice sites. The total-electron number is N = N a and the spin density is given by 
m = n-f — ni- 

The critical exponents which determine the response-function properties and the associate instabilities of the many- 
electron problem (1) are fully controlled by the pseudoparticle interactions. Let us then provide some basic information 
on the low-energy pseudoparticle operator basis which diagonalizes the quantum problem (1). This diagonalization 
uses the BA H||. We consider all finite values of U, electron densities < n < 1, and spin densities < m < n. (In 
Sec. V we also consider the n = 1 half-filling and m = n fully-polarized ferromagnetic phases.) For this parameter 
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space the low-energy physics is dominated by the lowest-weight states (LWS's) of the spin and 77-spin algebras which 
in Refs. [0-^9) were called of type I - in this paper we call them for LWS's I. These are described by real BA 
rapidities, whereas all or some of the BA rapidities which describe the LWS's II are complex and non-real Both 
the LWS's II and the non- LWS's out of the BA solution have energy gaps relative to each canonical ensemble ground 
state |r^-|l9| and do not contribute to the quantum problem instabilities studied in this paper. 

In the Hilbert subspace spanned by the LWS's I the BA solution was shown to refer to an operator algebra which 
involves two types of pseudoparticle creation (annihilation) operators b q a (b q . a ). These obey the usual anti-commuting 
algebra 

{t>l, a , bq', a '} = <W<W> { 6 g,a> & g',a'} = °> = • ( 2 ) 

Here a refers to the two pseudoparticle colors c and s 17|-19|. The discrete pseudomomentum takes on values 
qj = jf-If, where in contrast to the usual momentum, are consecutive integers or half-odd integers. It can be 
shown that in spite of the above anti-commuting algebra the a pseudoparticles have neither a fcrmionic nor a bosonic 
statistics [^T|. However, provided we introduce the topological-momentum shift operators studied in Refs. ||TT| , ^9| , pT| 
the second-quantization pseudoparticle representation can be used. There are dp a values of Jj* (dp a is the a-band 
Fock-space dimension), i.e. j = 1, ...,d,Fa- (The a topological-momentum shift operators produce shifts of in 
the a-band pseudomomentum numbers I" [ jTl||l9|]2l|| .) A LWS I is specified by the distribution of N a occupied values, 
which we call a pseudoparticles, over the dp a available values. There are dp Q — N a corresponding empty values, 
which we call a pseudoholes These are good quantum numbers such that dp c = N a , N c = + dp s = iV-j- 
and N s = N±. 

The numbers I? are integers (or half-odd integers) for N s even (or odd), and H are integers (or half-odd integers) 
for dp s odd (or even) [0. Therefore, all excitations which change N s — N± and dp s — N-\ by an odd integer 
always involve topological-momentum shifts |l^Jl^,|2l|. All the LWS's I can be generated by acting onto the vacuum 
\V) (zero-electron density) suitable combinations of pseudoparticle operators 0,^8) and are Slatter determinants of 
pseudoparticle levels of the simple form 

\cb;N C: N s )= J] {Y[blJ\V). (3) 

a— c,s q 

We emphasize that all LWS's I are characterized by fixed values of N s and dp s and therefore the pseudomomentum 
takes on numbers I" of well defined integer or half-odd integer character. The ground state (GS) of a canonical 
ensemble with N c and N s pseudoparticle numbers has the particular form 

q M 

\GS;N C1 N S )= J] [ f{ b\j\V). (4) 

a=C ' S 3=4"J 

When N a is odd (even) and the numbers I" are integers (half-odd integers) the pseudo-Fermi points are symmetric 
and given by 0@ <Tf1 = -<If1 = iFa - ^ where 

qFa = ^u- (5) 

On the other hand, when N a is odd (even) and Jj* are half-odd integers (integers) we have that — QFa an d 
— QfI = 1 Fa — W~ or 1f2 = 1 Fa ~~ W~ an< ^ ~<lFa = iFa- Similar expressions are obtained for the pseudo-Brioullin 
zones limits if we replace in the above expressions N c by N a . On the other hand, the pseudo-Brioullin zones 
limits q^' are always symmetric and given by qi + ^ = —qi — -^[dp s — 1]. 

In the pseudoparticle basis spanned by the LWS's I and in normal order relatively to the GS (4) the Hamiltonian 
(1) has the following form 0{l§ 



:H:=Y^HW, (6) 

where, to second pseudoparticle scattering order 



i=l 
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q . a 

^ = ^EE \f™>M) ■ #»(?) ■■■ N a ,(q>) : . (7) 

q,a q',a' 

Here (7) are the Hamiltonian terms which are relevant at low energy ]is[j . Furthermore, at low energy and small 
momentum the only relevant term is the non-interacting term . Therefore, the c and s pseudoparticles are non- 
interacting at the small-momentum and low-energy fixed point and the spectrum is described in terms of the bands 
e a (q) studied in Ref. Jl3[ . 

At higher energies and (or ) large momenta the pseudoparticles start to interact via zero-momentum transfer 
forward-scattering processes of the Hamiltonian (6) — (7). As in a Fermi liquid, these are associated with / functions 
0,0 which read 

faa'(q,q') = 2Kv a (q)$ aa ,(q,q') + 2Trv a ,(q')$ a , a (q' ,q) 

+ ^ E 27rV a"^a"a{jqFa",q)^ a "a'(jqFa",q') , (8) 
j=zkl a" — c,s 

where the pseudoparticle group velocities are given by v a (q) = an d v a = v a {qF a ) are the pseudo- Fermi group 

velocities. In expression (8) $ aa ' (q, q') mesures the phase shift of the a' pseudoparticle of pseudomomentum q' due 
to the zero-momentum forward-scattering collision with the a pseudoparticle of pseudomomentum q. These phase 
shifts describe the pseudoparticle interactions and are defined in Ref. (TJ]. They control the low-energy physics. For 
instance, the related parameters 

e aa . = <W + <*w (q { F« , ?F + i ) + (- iqpl , q ( il> ) , 3 = 0A, (9) 

play a determining role in the expressions we evaluate in this paper. There is a simple and direct relation between 
the pseudoparticle-interaction parameters and the quantities of CFT. For instance, the anti-symmetric combinations 
of two-pseudoparticle phase shifts (9), are nothing but the entries of the transpose of the dressed-charge matrix 
||. Both the symmetric [j = in Eq. (9)] and anti-symmetric [j = 1 in Eq. (9)] combinations of two-pseudoparticle 
phase shifts (9) fully control the instabilities of the quantum problem, as we find in the next sections. 



III. THE INSTABILITIES OF THE QUANTUM PROBLEM 

In order to study the divergences of the low-frequency response functions Rex 15 (k, oj) we start by characterizing the 
relevant GS transitions which determine these divergences. 

Following Ref. |Q , at low energy there are independent right iV+ and left N~ a pseudoparticle conservation laws in 
the Hubbard chain. Here iV+ (or N~) are the number of a pseudoparticles, with q > (or q < 0). While fixed values 
of N c and N s (and then of iVj and N± electron numbers) define a canonical ensemble, we introduce the concept of sub 
canonical ensemble which refers to fixed values of iVjr, N~, , and N~ . Obviously, the same canonical ensemble is 
realized by several sub-canonical ensembles. 

Since all LWS's I of form (3) are characterized by fixed values of the pseudoparticle numbers iV+, N~ , N+, and 
N~, we can associate with each sub-canonical ensemble a low-energy Hilbert subspace spanned by these LWS's I 
which have the same pseudoparticle numbers. We call pseudo-ground state (PGS) the LWS I (or LWS's I) of minimal 
energy in each of these Hilbert sub spaces. All PGS's are of the form 

\PGS;D C ,D S )= [] I II blJ\V), (10) 

where 

D a - AiV+ = -Ajv- (11) 

refers to the AAfJ changes relative to the GS (4) with the same N c and _/V s numbers. 

An important point is that the transitions between a GS and any LWS I can be separated into two types of 
excitations: (a) a GS - GS or GS - PGS transition which involves changes in the pseudoparticle numbers and (b) a 
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Landau-liquid excitation associated with pseudoparticle-pseudohole processes relative to the final GS or PGS which 
do not change the right and left pseudoparticle numbers. 

A remarkable property is that the low-frequency correlation-function divergences are fully determined by the above 
transitions (a), followed by transitions (b). This is because these divergences are determined by the quantum overlap 
between GS's and PGS's and the correlation-function excitations. The momentum of these transitions (a) is given by 



k = *^2 D ^Qi 



(12) 



where for the pure GS - PGS transitions (by pure GS - PGS transitions we mean here those which do not change the 
N c and N s pseudoparticle numbers) D a — 0,±1,±2, ... are always integers, whereas for the GS - GS transitions D a 
can either be integers or half-odd integers. 

The pseudoparticle-pseudohole processes (b) increase the value of the critical exponents which control the response- 
function divergences by positive integers. This leads to positive exponents associated with irrelevant correlation- 
function terms. Therefore, the divergences occur always precisely at momentum values fc = fco (see Eq. (12)) which 
correspond to GS - GS or GS - PGS transitions (a). 

In order to find the divergences of all low-frequency response functions Rex (fc, u) we have followed the steps of 
Refs. Jl(],[ll]] and combined our suitable generator pseudoparticle analysis [|l7 18 with the results obtained from CFT 
[0,0. The asymptotic expression of the correlation function % , '(x,t) in x and t space is given by the summation of 
many terms of form (3.13) of Ref. with dimensions of the fields suitable to that function. For small energy the 
corresponding correlation functions in k and u> space are obtained by the summation of the Fourier transforms of 
these terms, which are of the form given by Eq. (5.2) of Ref. M. However, the results of Refs. MM do not provide the 
specific expression at k = ko and small positive to. In this case the above summation is equivalent to a summation in 
the suitable final GS's or PGS's which correspond to different values for the dimensions of the fields. 

As was already referred in Sec. I, we emphasize that expression (5.7) of Ref. @] is not valid in our case. While 
we need the general response- function oj dependence at k — k , expression (5.7) of Ref. (?J is only valid for a very 
particular relation between (fc — fco) and u>. Moreover, the limit (fc — fco) — ► of that expression does not provide the 
general low-w behavior at fc = fco. 

Following Ref. @, we have solved the following general integral 



g(k ,ui) = dx dte 1 ' 

JO J -oo 



n 



(x + IV a t) 



2A< 



(13) 



which is the integral of expression (5.2) of Ref. at fc = fco. In Eq. (13) with t — +1 and b = — 1 for right and 
left a pseudoparticles, respectively, are the GS - GS or GS - PGS suitable dimensions of the fields §,[l7| of the general 
form 



2Ai 



AN a , 



1 2 



(14) 



where £^ Q , (with j — 0, 1) are the combinations of two-pseudoparticle phase shifts (9) and AN a refers to the changes 
in the numbers N a associated with the transition. We find for the integral (13), g(ko,ui) oc uo^ . Here 



so 



= 2^A^-2. 



(15) 



Comparing our expression with expression (5.7) of Ref. [j7| we confirm these expressions are different. This is because 
these expressions correspond to different limiting cases which do not commute. We find that at the momenta fc — fco 
[given by Eq. (12)] and low frequency u> all correlation functions behave as Imx^fco, u>) oc w w with the corresponding 
real part given by 



Rex (fco, oc cj ? 



(16) 



for ^ and by Rex (fco, w) oc — ln(w) for ^ = 0, respectively. Given this related expressions for the imaginary 
and real parts, below we consider only the latter part. Our task is detecting all response functions and momenta 
values fc = fco for which the corresponding exponent Eq. (15), is zero or negative. We have evaluated it for all 
correlation functions and possible values of fco. All divergences occur only for some of the correlation functions which 
show logarithmic divergences at U = 0, ie all divergences occurring at finite values of U correspond to logarithmic 
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divergences at U = 0. However, we find that the onsite electron-electron interactions supress such divergences 
in the case of the the superconductivity functions and for some values of U in the case of the charge and spin 
functions Rex p (±2fcpf, uS) and Rex a * (±2fcp j, uj), respectively. In general the effects of electron correlations replace 
the logarithmic divergence by a Luttinger-liquid [^of power-law divergence. 

The use in Eqs. (15) and (16) of suitable dimensions of the fields (14) for the different correlation functions leads 
to the following results. For the transverse-spin function we find 

Rcx^ {±2k F , uj) cx lu^ , (17) 

where the exponent is such that — 1 < <r s j_ < and given by 

?s± = -2 + 2^[(^) 2 + (^) 2 ]. (18) 

a 

The divergences of the charge and spin functions were already studied in Ref. |lC[ ] and are controlled by the same 
exponents. The expressions read 

Rex P (±2fc F(T , uj) cx w 5 ™ , Rex^ (±2k Fa , uj) cx w 5 ™ , (19) 

where the exponents are such that —\< <i C si < and — h < <r cs | < 2 and are given by 

w = -2 + 2^(cL) 2 , w = -2 + 2]T(£L-£L) 2 - (20) 

a a 

(We emphasize that at the parameter-space points which correspond to c cs f = the expressions of Eq. (19) are not 
valid and should be replaced by logarithmic uj dependences.) 
In the case of the singlet-superconductivity function we find 

Re X ss (±2k F ,Lu) cx w f "+ , Rex ss (±[fc^ T - k Fi ], uj) cx uf— , (21) 

where the exponents are such that < s aa + < 1 and < ? ss _ < 3 and are given by 

W = -2 + 2 £[(^) 2 + (C + %) 2 ] 



= -2 + 2 ]T[(^ e as ) 2 + (Cc + %) 2 ] ■ (22) 

a 

Since for U > these exponents are positive (as the figures of next section confirm) it follows that the electron- 
electron interactions remove the U = singlet-superconductivity instability. The same occurs with the triplet a 
superconductivity function which reads 

Rex* SCT (0,w)ocw ft <", (23) 
where the exponent is such that < q s | < 2 and < Cts| < 1 &ud is given by 

, tsl = -2 + 2 + Cs) 2 , ftst = -2 + 2 J2(£c) 2 ■ (24) 

a a 

The a one-electron Green function was already studied in Ref. jll| and for small uj and U > is such that 

ReG a (±k Fa ,uj) cx w fa , (25) 
where the exponent — 1 < ? CT < — 1/2 is given by 

* = - 2 + E ^(£c - £.) 8 + (e° c ) 2 ] , q = -2 + £ i[(£L) 2 + (£°c + e° s ) 2 ] ■ (26) 

a a 

All exponents (18), (20), (22), (24), and (26) have a Luttinger-liquid non-classical character being U, n, and m 
dependent functions. They are fully controlled by the two-pseudoparticle forward-scattering collisions [see Eqs. (9) 
and (16)]. 
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IV. COMPARATIVE STUDY OF THE INSTABILITIES 



In this section we present and discuss the U, electronic density, and magnetic-field dependence of the exponents 
(18), (20), (22), (24), and (26) which control the quantum-problem instabilities. Our comparative study of these 
exponents allows finding the dominant instabilities in the different regions of parameter space. 

In the Table we present various limiting values of the correlation-function exponents. There the m — > n values are 
also valid when H — > H c , U — > U c , and n — > n c , where H c , U e , and n c are the (interrelated) critical values for the 
onset of "ferromagnetism" (i.e., full spin polarization, m — ► ri) and are defined by Eq. (2) of Ref. Jl6| . Below we 
often use the normalized magnetic field h — H/H c . In the limit n — * 1 (not presented in the Table) the exponent 
expressions also simplify because = £° c = 1 and = £cs = and the remaining parameters change from = 1/2, 
e ss = 1/V2, £ s ° c = -1/^2, and e ss = V2 when m -> to & = 0, £, = 1, & = 0, and & = ; 1 as m - n. 

At finite values of U the exponents c s _i_, and <j cs | are always negative, whereas <r cs f is only negative for some 

regions of parameter space. This exponent is negative in the limit of m — > for all values of the density and U. 
For < m < n it can have both negative and positive values |l0[ . It is an increasing function of the magnetic field, 
becoming positive for a value of that field which depends on n and U. In the limit m — ► n it becomes positive for 
all finite values of U and all electronic densities n. For finite values of U all the remaining exponents are positive. 
In the figures we use often the electronic density n — 0.5 because it is typical for many of the quasi-one-dimensional 
Beckgaard salts j^g. 

In Fig. 1 the exponents (a) (b) fej and C C sT> ( c ) anc ^ (°0 ?T ano - m Fig- ^ the superconductivity exponents 
(a) 5ss+, (b) <r ss _, (c) q s |, and (d) ft s f are plotted vs U for n = 0.5 and various values the of the magnetic field. 

While the exponents s s ± and q, s ^ are always a decreasing function of U and the exponent c C s| a deacreasing and 
increasing function of U for smaller and large values of h, respectively, all remaining above exponents are always an 
increasing function of U. 

In Fig. 3 the exponents (a) <; s ^l, (b) <; cs | and ? C sT> ( c ) an d (d) ?f and in Fig. 4 the superconductivity exponents 
(a) ?ss+) (b) ? ss — , (c) q s |, and (d) ?t s | are plotted vs the electronic density for U — 10 and three values of the magnetic 
field. 

All the exponents are smooth functions of the density showing either a maximum or a minimum for an inter- 
mediate value of the density. While the exponents c s _L and <^ csa show a maximum value, the Green-function and 
superconductivity exponents show a minimum value. 

The magnetic-field dependences of the correlation-function exponents are illustrated in Fig. 5 and 6 where these 
parameters are plotted as functions of h for density n — 0.5 and various values of U. While the exponents s s ± and 
Stsi are always a decreasing function of the magnetic field, all the remaining exponents plotted in the figures are in 
general an increasing function of h. 

We emphasize that the exponent ?j is always smaller than or equal to the exponent In order to understand 
the parameter-space critical line defined by the equation <; s _l = ?| , the exponents Csi. and ^ associated with the two 
dominant instabilities are plotted in Figs. 7 (a)-(d) as functions of the magnetic field h and for several values of the 
electronic density n and onsite interaction (a) U = 1, (b) U = 3, (c) U = 5, and (d) U = 20. At zero magnetic field 
?! = sj, < <; s _\_, and the exponent <; s _l equals the exponents <; cs ± and £ cs f which are equal in that limit. (The latter 
were already studied in Ref. [fl0|| .) Therefore, adding or removing of electrons is in this case the dominant instability. 
(In the case of the weakly coupled-chain Hubbard system this leads to a dominance for interchain electron hopping.) 
Also at small values of U that instability remains dominant for all values of the electronic density and magnetic field. 

However, the magnetic field removes the SU (2) spin symmetry and for intermediate and large U values induces a 
dominant character for the spin-flip excitation over adding or removing of single electrons. This occurs for U > U*, 
where U* — U*(n,h) is a critical value. Depending on the parameters that we keep constant, we can also define 
critical values n*(U,h) and h*(n, U). The critical values are obtained by solution of the above equation, = Cs-Lj 
upon suitable boundary conditions. Therefore, they correspond to the points of Figs. 7 (a)-(d) where the two exponent 
lines cross. As these figures confirm, U* changes from 

4<sin(7m) 

U* = )=-!■ — , 27 

tan(§[V2-l]) 

when h — ► 1 to U* = oo when h — * 0. 

In Figs. 8 and 9 the functions U*{n) and h*(n) are presented for constant values of h and U, respectively. These 
define the critical lines which divide the regions of parameter space where |-electron adding and removing on the one 
hand and spin-transverse SDW on the other hand are the dominant instability. 

At finite values of the magnetic field h < 1, |-electron adding or removing remains dominant for U < U* (see Fig. 
8). This region is larger for intermediate densities and for small values of the magnetic field and smaller for small 



7 



densities and densities close to half filling and for large values of h. Therefore, the function U*(n) plotted in Fig. 8 
shows a maximum at intermediate electronic densities and is larger at smaller values of the magnetic field h. 

On the other hand, at intermediate and large values of U, "["-electron adding or removing remains dominant for 
h < h* (see Fig. 9). This region is larger for intermediate densities and for small values of U and smaller for small 
densities and densities close to n = 1 and for large values of U. It follows that the function h*(n) plotted in Fig. 9 
shows a maximum for intermediate values of the electronic density and is larger for smaller values of U. 



V. THE METAL - INSULATOR HALF-FILLING AND FULLY-POLARIZED FERROMAGNETIC 

TRANSITIONS 

In this section we consider only these correlation functions whose exponents are negative for finite U and omit the 
study of the superconductivity correlation functions. 

When we change the chemical potential [i from fi > A ^ H to fj, < A " H , where A mh is the half- filling Mott-Hubbard 
gap H which at finite magnetic field was studied in Ref. Jl3| (and thus the electronic density changes from n < 1 
to n — 1), there occurs in the many-electron system a metal - insulator transition |16[ ]. Also when we change the 
magnetic field from H < H c to H > H c (and thus the spin density m changes from m < n to m = n), there occurs 
a transition to fully-polarized ferromagnetism. In this section we find that i)- some exponents show a discontinuity 
because they have different values in both sides of the above transitions and ii)- these transitions remove some of the 
low-energy power-law divergences in the correlation functions but in some cases these divergences are shifted to finite 
energies. 

At half filling both the charge and the one-particle excitations associated with creation of one electron cost a 
minimal energy Amh- Therefore, both the charge response function and the one-electron Green function vanish for 
positive energy, oj < Amh- 

On the other hand, there are gapless spin excitations which lead to low-energy power-law behavior for the spin 
response functions. However, the half-filling metal - insulator transition leads to discontinuities in the corresponding 
exponents. We find 

Rex^ (±tt, ui) oc w Ssi , Rex az {±2k Fai uj) cx uf"" , (28) 

where the exponents c s j_, Scsl, ancl CcsT are f° r n —* 1 and at n = 1 (and for // — > Aa 2 fH from above and for \jl < A, ^ H ) 
given by 

s s ± = -2[§ - (-^r) 2 ] , ^4 = -2[i - (e cs ) 2 - (d) 2 ] , <-t - -2[i - (i - e cs r - (e ss ) 2 } , (29) 



and 

S S ± = -2[1 - (^r) 2 ] , W = -2[1 - (CL) 2 ] , fe T = -2[1 - (CL) 2 ] , (30) 

respectively. In the limit m — > these exponents read 

?s-L = fesj, = ^csj = — 1/2 ! (31) 

for n — * 1 (and for /i — > A ^ rff from above) and 

Ss± = Tcsl = <TcsT = — 1 ' (32) 

for n = 1 (and for /i < A f ii ). 

On the other hand, in the limit m — > n these exponents are given by 

T s ± = -1 , fcsj. = , Test = 2 , (33) 

for n — > 1 (and for /i — > Amjl from above) and 

? s± = -3/2 , W = , ? csT = , (34) 
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for n — 1 (and for /i < ^ H ). Obviously, in the cases where the exponents vanish expressions (28) should be replaced 
by a logarithmic function. Note that the logarithmic behavior is only reached in the limit m — ► n since otherwise the 
correspondent exponents are different from zero. 

In the case of fully-polarized ferromagnctism there are no spin-down electrons in the system. In this case a spin 
flip costs a minimal energy given by 

A F = 2 tM) [H-H c }, (35) 

where H > H c and A F is the fully-polarized ferromagnetism gap. Therefore, the transverse-spin response function 
and the down-spin one-electron Green function vanish for positive energy, to < A F . Also, there is no k = kp[ = 
divergence in Kex p (k,uj) and Rex a *(k,uj). 

The charge response function, spin response function Rex CTz , and up-spin one-electron Green function all have a 
low-energy power-law structure in both sides of the ferromagnetic transition. That transition leads to discontinuities 
in the exponents of these correlation functions, which read 

Rex P (±2fc F ,w) oc w ScsT , Rex az {±2k F ,Lj) cx ui^'t , ReG T (±2fc F , u>) oc cu^ , (36) 

where for m —* n (and H — > H c from below) the exponents are given by 

e csT = 2[1 - vo? , ^ T = -I + i[l-77 ] 2 , (37) 



and 

'At sin(7rn) 



/ 2 n -if 1 
?7o = (-)tan I- 



U 



(38) 



On the other hand, atm = n (and H > H c ) the above charge and spin response functions have a non-interacting 
logarithmic behavior 

Rex P (±2fc F ,w) cx -ln(w), Re% CT * (±2fc F , u) oc - ln(w) , (39) 

because 



fe T = , (40) 

and the up-spin Green function exponent reads 

? T = -1. (41) 

Based on the occurrence of only zero-momentum pseudoparticle forward scattering and on the pseudoparticle- 
number conservation laws pH we find that in the case of half filling there are power-law divergences at finite energies, 
for instance, in the one-electron Green function. Finite-energy power-law divergences also occur, for example, in the 
transverse-spin response function in the fully-polarized ferromagnetic phase. (The above property of the pseudoparticle 
interactions is valid at all energy scales |2l| . A complete study of the quantum problem for all energies involves an 
infinite number of pseudoparticle branches - fortunately, the finite-energy divergences we study in this paper only 
involve the c and s branches.) 

By generalizing the pseudoparticle operator basis to half filling and the fully-polarized phase, we find an equivalent 
low-energy pseudoparticle theory where some of the low-energy u) divergences of Sec. Ill are replaced by low-energy 
(uj — Amh) and (w — A F ), respectively, divergences. At these energies the half-filling and fully-polarized ferromagnetic 
pseudoparticle theories involve the same parameters as the general < n < 1 problem studied in previous sections. 
For instance, for (u) — Amh) > and (oj — A F ) > all pseudoparticle-Hamiltonian parameters are obtained by 
considering the limits n-t 1 and m — ► n, respectively, of the corresponding < n < 1 and < m < n expressions 
(and not the corresponding n — 1 and m = n expressions). 

In the case of half filling the lack of charge gapless excitations and the associate opening of the Mott-Hubbard gap 
leads to the following expression for the one-electron Green function 

RcGcr (ifc Fcr , lu) oc [ui - A MH p , (42) 

where the exponent <; a is given by 
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?T =-2+i[(l- d) 2 + 1 + (£]J 2 + (d) 2 ] , (43) 

q = -2 + i[(d) 2 + 1 + (CD 2 + (d + O 2 ] , (44) 

and (w — Amjj) is small and positive. 

Both in the limits m — > and m — > n the exponents (43) and (44) for up-spin and down-spin electrons, respectively, 
are equal and read 

q=q = -7/8, (45) 

and 

q=?4 = -l/2, (46) 

respectively. 

On the other hand, in the fully-polarized ferromagnetic phase there are no gapless excitations in the transverse-spin 
channel. We find the following expression 

Re X a± {±2k Fl uj) oc [w- Af] 5 ' 1 , (47) 

where the exponent is given by 

< s ± = -I + M 2 ■ (48) 

Obviously, this expression refers to small positive values of the energy (w-Ap) and is not valid in the limit of U — > 
where <r s j_ — > 0. 

Equations (42) — (48) reveal that the opening of the Mott-Hubbard and ferromagnetic gaps associated with the 
electronic-correlation and magnetic effects, respectively, moves some of the uj — ► instabilities to lu — > Amh and 
u> — > Ap, respectively. These divergences occur at the same momenta values than the corresponding low- frequency 
< n < 1 and < m < n divergences. 

The finite-frequency expressions (42) — (48) cannot be derived by CFT. Their validity follows from the purely non 
dissipative character of the pseudoparticle interactions which survives at all energy scales of the quantum problem 



VI. COUPLED CHAINS 

The exponents calculated in the preceeding sections give the scaling dimension of the corresponding operators, in 
the presence of the interaction. Once they are known, we can analyze the effects of perturbations added to (Q) which 
can be written as a sum of terms containing the operators analyzed above. For simplicity, we first examine the case 
in which these operators couple two different chains. Let us write % = (O^ (x,t)O i} ). The Fourier thansform of 

X& goes at low frequencies as aA* . An interchain coupling term which involves the operator C 3 can be written as 

Wi,a = J# J dxOf(x)Of(x) , (49) 

where the labels 1 and 2 refer to the two different chains. We now add this term to (0) and study the scaling of as 
the cutoff is reduced. We define the dimensionless coupling J§ — J^/lu c , where uj c is an upper cutoff below for which 
the approximation of xn by a power law is approximately valid. Then, neglecting the influence of the other interchain 
couplings we find ]24|-|26| 

f) f 

lu c -^ = (-1 - 2») J,? . (50) 

OUJc 

Thus, the relevance of 0$ depends on whether q# is greater or less than —1/2. When < —1/2, 0$ is a relevant 
operator which grows upon scaling. The Renormalization-Group scheme used to derive ( fl9[ ) breaks down when J$ *~ 1. 
This defines a crossover scale, 
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Below this crossover the coupling J$ dominates the physics of the coupled-chain system and a new method needs to 
be applied. In the case of many-coupled chains, forming a 2D or 3D array, mean-field approaches are, most likely, 
adequate. Note that this analysis can be extended, in a straightforward way, to a many-chain system. This crossover 
has been studied for a variety of interactions in Ref. |27|] . 

The most commonly studied perturbation has been the interchain hopping. Then the required exponents are 
and which are given in Eq. (26). These operators are always relevant except as H — > H c . The fixed point which 
describes an isolated Luttinger liquid is destabilized and interchain phase coherence develops. Note, however, that as 
in the presence of the magnetic field ^ <?j,, the interchain hopping becomes spin dependent at low energies. This 
effect should influence the transport properties. It would be interesting to check this fact experimentally. 

So far, he have ignored the influence of the possible couplings among themselves. To leading order in the strength of 
the couplings, this is, at least, a second-order process, while the effect which leads to (|m]) is first order. It may happen, 
however, that new couplings are generated, which were zero at the initial stages of the renormalization process, i. 
e., in the bare Hamiltonian. When these new couplings are more relevant than the ones initially present, they may 
change drastically the physics at low energies. 

The coupling most likely to give rise to new interactions is the interchain hopping. It involves two electronic 
operators, and, to next order, it can generate four fermion operators, like interchain magnetic exchange, or Cooper 
pair hopping. In the following, we analyze the process by which an exchange term is generated if, initially, interchain 
hopping is present. We consider the transverse-spin exchange coupling, as this is the other most relevant operator in 
the presence of electron-electron repulsion, besides the hopping. From the calculations presented in the preceeding 
sections, it may become more relevant than the hopping itself in the presence of a magnetic field. The full scaling 
equation for this coupling, J s ±, is 

dJ 

Uc-^ = (-l-2>u±)J a ±-TJ t J l , (52) 

OUJ c 

where T is a constant of order unity. The second term in the right-hand side stands for the fact that a simultaneous 
hopping of an up-spin electron in one direction and a down-spin electron in the opposite generates a spin-exchange 
process. 

We can integrate Eqs. ( j50|j5lj ) with the initial conditions: Jf = J[ = Jo and J s ± — 0. The scale-dependent coupling 
J s j_ depends on / = log(w/w c ) as 

Js±(l) « —, L . (53) 

The low-energy behavior of the coupled system shows two different regimes, depending on the value of the denom- 
inator in the right-hand side term in Eq. (p3[): 

i) If the denominator is of order unity, the interchain hopping diverges first, at a scale given by Eq. (|5l|). Below this 
scale, a system of coupled chains resembles a 2D or 3D anisotropic Fermi liquid, with a residual exchange interaction. 
A mean-field analysis of this model leads to the possibility of a SDW ground state. 

ii) When the denominator is much less than one, the induced transverse-spin exchange grows faster than the 
hopping term. Then, before coherence between the single-electron wavefunction at the different chains is established, 
the exchange-term becomes of the order of the cutoff. The system probably develops a spin gap, as in the previous 
case, but the transport properties should be drastically different. 

Note, finally, that this analysis is valid for arbitrary strength of the repulsion, U, and weak interchain couplings. 
It would be interesting to connect this picture to the behavior obtained in the opposite limit of strong-interchain 
hopping and weak-interchain coupling, studied recently p8| , ^9| . 



VII. CONCLUDING REMARKS 



In this paper we have presented a comparative study of the Hubbard chain instabilities. In contrast to the one- 
electron ID quantum problem, the many-electron instabilities are not associated in general with logarithmic diver- 
gences in the low- frequency response functions. Instead these functions show Luttinger- liquid power laws which are 
controlled by non-classical exponents. Moreover, the singlet and triplet superconductivity instabilities are removed 
by the effects of the electron correlations. 
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We have studied in detail the U, electronic density n, and magnetic-field H dependence of all the exponents which 
are associated with instabilities of the many-electron quantum problem. The interplay of the magnetic field and 
electron correlations renders the transverse-spin SDW instability at momenta ±2fcp dominant for large regions of 
parameter space. By removing the SU (2) spin symmetry the magnetic field also makes the exponents associated with 
the spin and transverse-spin response functions different. 

We went beyond CFT and characterized divergences which occur at finite energy values in some correlation functions 
in the cases of the half-filling and fully-polarized ferromagnetic phases. We have also studied the discontinuities which 
occur in some power-law exponents which have different values in both sides of the transitions to or from the above 
phases. 

Our exact results might be useful for the study of the electronic instabilities in real quasi-one dimensional materials. 
For instance, transverse SDW's phases were often observed in real low-dimensional materials (2^j2^]. Following the 
discussion of Sec. VI, elsewhere we will use the present instabilities of the ID problem to construct a suitable effective 
Hamiltonian for the description of a system of weakly coupled Hubbard chains. 
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TABLE - Limiting values of the exponents whose expressions are defined in Eqs. (18), (20), (22), (24), and (26). 
Here the m — > parameter £o = £c C changes from £o = V2 at {/ = to £o = 1 as U — > oo and the parameter 770 is 
given in Eq. (38). 
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FIG. 1. The exponents (a) <; s ±, (b) ? cs ^ and ? cs |, (c) and (d) ?| defined by Eqs. (18), (20), and (26) as functions of U for 
n = 0.5 and for values of the magnetic field ft — 0.3 (full line), ft = 0.6 (dashed line), and ft = 0.9 (dashed-dotted line). 

FIG. 2. The exponents (a) <j S s+, (b) ? ss -, (c) Stsi, an d (d) CtsT defined by Eqs. (22) and (24) as functions of U for n = 0.5 
and for values of the magnetic field ft — 0.3 (full line), ft = 0.6 (dashed line), and ft = 0.9 (dashed-dotted line). 

FIG. 3. The exponents (a) q a ±, (b) ? cs | and Ccsf, (c) and (d) ?f defined by Eqs. (18), (20), and (26) as functions of the 
electronic density n for C7 = 10 and for values of the magnetic field h = 0.1 (full line), ft = 0.2 (dashed line), and ft = 0.3 
(dashed-dotted line). 

FIG. 4. The exponents (a) ? ss +, (b) ? ss _, (c) ? ts |, and (d) c tsT defined by Eqs. (22) and (24) as functions of the electronic 
density n for U = 10 and for values of the magnetic field h = 0.1 (full line), /i = 0.2 (dashed line), and ft = 0.3 (dashed-dotted 
line) . 

FIG. 5. The exponents (a) <; s ±, (b) ? cs | and ? cs |, (c) and (d) ?| defined by Eqs. (18), (20), and (26) as functions of the 
magnetic field ft for n = 0.5 and various values of U. 

FIG. 6. The exponents (a) (b) ? s *-, (c) ?t s |, and (d) ? ts f defined by Eqs. (22) and (24) as functions of the magnetic 

field ft for n = 0.5 and various values of U. 

FIG. 7. The exponents <; s ± and ?| associated with the two dominant instabilities as functions of the magnetic field ft and for 
several values of the electronic density n for (a) U = 1, (b) U = 3, (c) [7 = 5, and (d) U = 20. Comparison of the values of 
both exponents shows that for intermediate and large values of U the transverse SDW becomes dominant for magnetic fields 
ft larger than about 0.75 to 0.8. 

FIG. 8. The critical line for U*(n) in the U — n plane for values of the magnetic field ft = 0.8 (full line) and ft = 0.9 (dashed 
line) . 

FIG. 9. The critical line for ft*(n) in the ft — n plane for values U = 10 (full line) and U = 20 (dashed line). 
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